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Abstract 

Starting from the Liouville equation, we derive the exact hierarchy of equations sat- 
isfied by the reduced distribution functions of the single species point vortex gas in two 
dimensions. Considering an expansion of the solutions in powers of 1/N (where N is the 
number of vortices) in a proper thermodynamic limit N — ► +oo, and neglecting some col- 
lective effects, we derive a kinetic equation satisfied by the smooth vorticity field which is 
valid at order 0(1/N). This equation was obtained previously [P.H. Chavanis, Phys. Rev. 
E, 64, 026309 (2001)] from a more abstract projection operator formalism. If we consider 
axisymmetric flows and make a markovian approximation, we obtain a simpler kinetic 
equation which can be studied in great detail. We discuss the properties of these kinetic 
equations in regard to the i7-theorem and the convergence (or not) towards the statistical 
equilibrium state. We also study the growth of correlations by explicitly calculating the 
time evolution of the two-body correlation function in the linear regime. In a second 
part of the paper, we consider the relaxation of a test vortex in a bath of field vortices 
and obtain the Fokker-Planck equation by directly calculating the second (diffusion) and 
first (drift) moments of the increment of position of the test vortex. A specificity of our 
approach is to obtain general equations, with a clear physical meaning, that are valid for 
flows that are not necessarily axisymmetric and that take into account non-Markovian 
effects. A limitations of our approach, however, is that it ignores collective effects. 
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1 Introduction 



Several authors have wondered whether fluid turbulence could be described in terms of sta- 
tistical mechanics [T]. Three dimensional turbulence has been attacked by different methods 
[21 [3j HI [5] inspired by statistical mechanics and kinetic theories. Some progress has also been 
made in the simpler case of two dimensional turbulence (see reviews in [6j [7J [8]). Two dimen- 
sional turbulence is not just academic but is relevant to describe geophysical and astrophysical 
flows. Two dimensional flows are characterized by the spontaneous formation of large-scale vor- 
tices that dominate the dynamics [HI [10]. The most famous example is Jupiter's great red spot, 
a huge vortex persisting for more than three centuries in a turbulent shear layer between two 
zonal jets in the southern hemisphere of the planet [II]. Other examples of this self-organization 
are the cyclones and anticyclones in the earth atmosphere, the jets in the oceans like the gulf 
stream or the intense jets on Jupiter [12J. As a first step to tackle the problem, it can be of 
interest to study the dynamics of a system of A point vortices on a plane [13]. Each vortex 
produces a velocity field that moves the other vortices in a self-consistent manner. The velocity 
created by a vortex decreases like 1/r which is similar to the Coulombian or Newtonian interac- 
tion in two dimensions. Therefore, the interaction between point vortices is long-range, like the 
interaction between stars in a galaxy or between electric charges in a plasma. Note, however, 
that point vortices produce a velocity while material particles produce a force (acceleration). 
Apart from this (important) difference, the point vortex gas has a Hamiltonian structure [2] 
and we can try to apply the methods of statistical mechanics and kinetic theory to that system. 
Therefore, the A-vortex problem [13] is of fundamental interest in statistical mechanics and 
kinetic theory. It provides a physical example of systems with long-range interactions, whose 
dynamics and thermodynamics are actively studied at present [T5j . 

The statistical mechanics of 2D point vortices was first considered by Onsager [IB] in a sem- 
inal paper. He showed that statistical equilibrium states with sufficiently large energies have 
negative temperatures. For such states, like-sign vortices have the tendency to group them- 
selves and form clusters. If the circulations of all the point vortices have the same sign, the 
equilibrium state is a large-scale vortex (supervortex) similar to vortices observed in geophysical 
and astrophysical flows. When the point vortices have positive and negative circulations, the 
equilibrium state is generically a dipole made of a cluster of positive vortices and a cluster of 
negative vortices. The pioneering work of Onsager was pursued by Joyce & Montgomery [17] 
and Lundgren & Pointin [18], using a mean field approximation. Using a combinatorial analy- 
sis, Joyce & Montgomery introduced an entropy for the point vortex gas which is similar to the 
Boltzmann entropy for material particles. The statistical equilibrium state (most probable) is 
obtained by maximizing this Boltzmann entropy while conserving all the constraints imposed by 
the dynamics (total number A of point vortices and energy E, as well as angular momentum L 
and impulse P for domains with a special symmetry). For point vortices with equal circulation 
7, the smooth vorticity field is given by the Boltzmann distribution uj(t) = Ae~^^ r \ where the 
potential is played by the stream function ip(r). Using uj = —Aip, the stream-function is then 
determined by the Boltzmann-Poisson equation. Lundgren & Pointin started from the exact 
equilibrium hierarchy of equations satisfied by the reduced distribution functions Pjivi-, ■ 
of the point vortex gas and, by neglecting all the correlations between point vortices, derived a 
differential equation determining the equilibrium distribution of the one-body distribution func- 
tion Pi(ri). Using the fact that uj(t) = A7Pi(r), the mean field equation derived by Lundgren 
& Pointin coincides with the Boltzmann-Poisson equation derived by Joyce &; Montgomery. In 
a mathematical work, Caglioti et al. [IH] showed rigorously that the mean field approximation 
is exact in a proper thermodynamic limit A — > +oo such that 7 ~ 1/A, E ~ 1, (3 ~ A and 
V ~ 1 (where V is the area of the domain). In that limit the A-body distribution at statistical 
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equilibrium is a product P/v(r 1; ...rjv) = Pi(r 1 )...P 1 (r A r) of N one-body distributions that are 
solution of the Boltzmann-Poisson equation. This statistical equilibrium state is expected to 
be achieved for t — > +00. We stress, however, that the statistical theory is based on the as- 
sumption that "at statistical equilibrium, all accessible microstates are equiprobable" . This is 
essentially a postulate, so there is no guarantee that the point vortex gas will reach a statistical 
equilibrium state of the form described above (based on the microcanonical distribution). In 
order to determine the timescale of the relaxation of the smooth vorticity field u(r,t), and in 
order to establish whether (or not) the system will truly relax towards Boltzmann statistical 
equilibrium, we must develop a kinetic theory of point vortices. 

A kinetic theory was developed by Dubin & O'Neil [20] in the case of a non neutral plasma 
confined by a strong magnetic field, a system isomorphic to the point vortex gas. They started 
from the Klimontovich equation and used a quasilinear approximation to determine the current 
of the smooth density due to discrete interactions between point vortices. They considered an 
axisymmetric evolution of the system and, in the course of their derivation, made a Markov 
approximation assuming that the two-body correlation function relaxes on a timescale that 
is much shorter than the timescale on which the smooth density field changes (this is the 
counterpart of the Bogoliubov hypothesis in plasma physics) . They obtained a closed expression 
of the current, see Eq. (11) of [20], taking into account "collective effects" between the particles. 
These collective effects are similar to those giving rise to the Debye shielding in plasma physics 
in the Lenard-Balescu approach [211 122] . In plasma physics, they take into account the fact 
that a charge is surrounded by a polarization cloud of opposite charges. In the case of point 
vortices, their physical interpretation and their consequence is more difficult to establish. 

A kinetic theory of point vortices was carried out independently by Chavanis [23], using an 
analogy with the kinetic theory developed for stellar systems. He started from the Liouville 
equation and used the projection operator formalism of Willis & Picard [24J to derive a kinetic 
equation for the smooth vorticity distribution u(r,t). By this method, he obtained a kinetic 
equation, see Eq. (128) of [23], that is valid for flows with arbitrary symmetry (non necessarily 
axisymmetric) and taking into account memory effects. This is the counterpart of the general- 
ized Landau equation in stellar dynamics derived by Kandrup [25] using the same formalism, 
see Eq. (42) of [25]. If we restrict ourselves to axisymmetric flows and make a Markovian 
approximation, this leads to a simplified kinetic equation, see Eq. (133) of [23], which coin- 
cides with the equation obtained by Dubin & O'Neil [2U] when collective terms are ignored. 
In a sense, the simplified kinetic equation (133) obtained by Chavanis [23J is the counterpart 
of the Landau [26J equation in plasma physics while the more general equation (11) obtained 
by Dubin & O'Neil [20] is the counterpart of the Lenard-Balescu equation [211 122]. Finally, 
the general non-Markovian equation (128) of [23] is related to the Master equations discussed 
by Prigogine [27] in plasma physics. Therefore, there are many interesting analogies between 
plasma physics, stellar dynamics and vortex dynamics. In these analogies, the position r of the 
point vortices plays the role of the velocity v of the particles in a plasma or in a stellar system, 
and the angular momentum L = J ur 2 dv plays the role of the kinetic energy K = | J fv 2 dv. 
Chavanis [23] also considered the evolution of a test vortex in a bath of field vortices with fixed 
distribution (e.g. a thermal bath at statistical equilibrium) and used the projection operator 
formalism to derive a Fokker-Planck equation for the evolution of the one-body distribution of 
the test vortex P(r, t) in the bath. This equation involves a term of diffusion and a term of 
drift that are both position dependent. For a thermal bath, i.e. when the field vortices are at 
statistical equilibrium, the diffusion coefficient and the drift coefficient (mobility) are related to 
each other by a sort of Einstein relation involving a negative temperature (in cases of physical 
interest). The resulting Fokker-Planck equation, see Eq. (115) of [23], is the counterpart of the 
Kramers-Chandrasekhar equation, see Eq. (10) of [28], in stellar dynamics. 
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The projection operator formalism which is at the basis of the above-mentioned kinetic 
theory is very abstract and it is not clear to understand which approximations have been made 
in the course of the derivation. In this paper, we present an alternative derivation of the 
kinetic equations obtained in [23] which is more transparent. We start from the exact out-of- 
equilibrium BBGKY-like hierarchy of equations satisfied by the reduced distribution functions 
Pj(ri, Tj, t) of the point vortex gas and consider an expansion of the solutions in powers 
of 1/N in a proper thermodynamic limit N — > +00. The kinetic equations obtained in [23] 
are recovered at the order 0(1/N). The idea of using a BBGKY-like hierarchy to derive the 
kinetic equations of [23] was given in [29]. This derivation has the advantage of being much 
simpler than the previous one and shows clearly the domain of validity of the kinetic equations. 
It also shows how collective terms can be included in the calculations. However, we shall not 
try to evaluate these terms in the present paper because they require heavy calculations in 
the complex plane similar to those performed in plasma physics to derive the Lenard-Balescu 
equation from the BBGKY hierarchy. We hope to come to this problem in a future work. 

The paper is organized as follow. In Sec. [2] we consider the statistical equilibrium state. 
We recall the equilibrium hierarchy derived by Lundgren & Pointin [IB] and the proper ther- 
modynamic limit for the point vortex gas. For N — > +00, the two-body correlation function 
vanishes so that the mean field approximation is exact in that limit. We derive the differential 
equation satisfied by the smooth vorticity field. For the usual potential of interaction between 
point vortices, it reduces to the Boltzmann-Poisson equation. We go beyond the mean field 
approximation and determine the exact integrodifferential equation satisfied by the two-body 
correlation function at order 0(1/N). In Sec. [31 we consider the out-of-equilibrium problem 
and derive the BBGKY-like hierarchy satisfied by the reduced distribution functions of the sin- 
gle species point vortex gas. We close the hierarchy by considering an expansion of the solutions 
in powers of 1/N in the proper thermodynamic limit N — > +00 (Sec. 13. ip . To leading order in 
N — » +00, we find that the smooth vorticity field satisfies the 2D Euler equation. This is the 
counterpart of the Vlasov equation in plasma physics and stellar dynamics. At order 0(1/N), 
we obtain an exact system of coupled equations taking into account "distant collisions" between 
point vortices. If we neglect collective effects, we can obtain an explicit kinetic equation that is 
valid for flows that are not necessarily axisymmetric and that takes into account non-markovian 
effects. It is valid on a timescale ~ Nto- We simplify this equation by considering axisymmetric 
flows and arguing that markovian effects can be neglected for N — > +00 (Sec. I3.2p . We discuss 
the properties of these kinetic equations in regard to the if -theorem and the convergence (or 
not) towards the statistical equilibrium state (Sec. 13. 5p . We also study the growth of corre- 
lations by explicitly calculating the time evolution of the two-body correlation function in the 
linear regime (Sec. 13.61) . In Sec. HJ we show that the above-mentioned kinetic equations can 
also be obtained from a quasilinear theory starting from the Klimontovich equation. This is 
the approach considered by Dubin & O'Neil [2D] for axisymmetric flows. We show how it can 
be generalized to arbitrary flows when collective effects are neglected. We also stress the con- 
nection with the quasilinear theory of the 2D Euler- Poisson system developed by Chavanis [30] 
to describe the process of violent relaxation in the collisionless regime [3T] [32] [33] 134"] . Finally, 
in Sec. \5[ we consider the relaxation of a test vortex in a bath of field vortices at equilibrium 
and obtain the Fokker-Planck equation by directly calculating the second (diffusion) and first 
(drift) moments of the increment of position of the test vortex. We obtain general expressions, 
with a clear physical meaning, that are valid for flows that are not necessarily axisymmetric and 
that take into account non-Markovian effects. We simplify them in the case of axisymmetric 
flows and recover by a direct calculation the Fokker-Planck equation that was obtained in [23] 
from the more formal projection operator formalism. We also indicate how the results can be 
generalized when collective effects are taken into account in the problem. 
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2 The statistical equilibrium state 



To set the notations and show the connection with the kinetic theory developed in the next 
section based on a BBGKY-like hierarchy, we first derive the differential equation satisfied by 
the smooth vorticity profile at statistical equilibrium from an equilibrium hierarchy [T8] . 



2.1 The equilibrium hierarchy 

The exact vorticity field of a gas of point vortices with equal circulation 7 is given by 

N 



= "*<(*)), (1) 



where rj(t) is the position of the i-th point vortex at time t. The dynamical evolution of the 
point vortices is governed by the Hamilton equations 

dxi dH dxji dH 

^ dt di/i ' ^ dt dxi ' 

i/ = 7 2 ^n(|r,,-r J |), (2) 

i<j 

where the positions (x, y) of the point vortices are canonically conjugate [Hj. For simplicity, 
we present the results in an unbounded domain where the potential of interaction depends only 
on the absolute distance between point vortices, but we stress that most of our results can be 
extended to bounded domains by using the generalized Green function of Lin [35] satisfying 
the reciprocity property u(rj, r_j) = u(r,-,rj). The N- vortex distribution function satisfies the 
Liouville equation 

i=l 

expressing the conservation of the probability density P/v(ri, tn, t) that the first point vortex 
is in ri, the second in r-2 etc., at time t and where V, is the total velocity of point vortex %. It 
is created by the other vortices, so that 

V< = X)V(j-0, ( 4 ) 

where 

V(j - i) = - 7 z x (5) 

is the velocity created by point vortex j on point vortex %. We shall essentially consider the 
standard potential of interaction Uy = —(1/2%) In |r»— Tj \ but we leave the function u(\Ti— Tj\) as 
general as possible in order to describe different models like, e.g., the quasi-geostrophic model. 
Any function of the constants of motion of the Hamiltonian dynamics (energy E = H, angular 
momentum L = 7 ^ rf if the domain has rotational symmetry and impulse P = 7 £\ rj if the 
domain has translational symmetries) is a stationary solution of Eq. ([3]). For brevity, we shall 
only consider the conservation of energy (the case of an infinite domain with conservation of 
angular momentum is treated in [18]). The basic postulate of statistical mechanics states that, 
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at equilibrium, all microscopic configurations that are accessible (i.e. that have the correct 
value of energy) are equiprobable. There is no guarantee that the dynamics will lead the system 
to that "uniform" state because we could imagine that some regions of the 2iV-dimensional 
phase space could be more probable than others. However, if we accept this postulate, the 
equilibrium iV-body distribution is given by the microcanonical distribution 



P N (ri, ...,r N ) 



9(E) 



S[E-H{n,...,T N )]. 



(6) 



Using the normalization condition J P^ \\ i dvi = 1, we deduce that the density of states with 
energy E is given by 

g{E)= I S[E-H(r u ...,r N )]l[dT i . (7) 

The microcanonical entropy of the system is defined by S(E) = \ng(E) and the microcanonical 
temperature by 1/T(E) = dS/dE (we take the Boltzmann constant ks = 1). We introduce 
the reduced probability distributions 



Pjv(ri, ...,r N )dr j+1 ...dr 



N- 



(8) 



For identical particles, the smooth (average) vorticity field is related to the one-body distribu- 
tion function by 



iV 
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» = <^ 7( J(r-r i )>=iV7P 1 (r). 



(9) 



i=i 



Note that the vorticity field is proportional to the density of point vortices: oo(r) = jn(r). The 
total circulation is T = J u(r)dr = and the average value of the energy is 

E = (H) = ^N(N - 1) 7 2 J u(\r- r'|)P 2 (r, r')drdr f . (10) 

), we can obtain an equilibrium 



By differentiating the defining relation for Pj and using Eq. 
hierarchy of equations for the reduced moments [18]: 







dPi 



dr x g(E) dE 



9(E)P t 



E 

i=2 



7 



{N-jW 



1 d 



dr 1 g(E) dE 



9{E)P j 



;+i 



dr 



j+i- 



'ir 



This is the counterpart of the equilibrium hierarchy in plasma physics. It is however more 
complex in the present situation because it has been derived in the microcanonical ensemble. 
Since statistical ensembles are generically inequivalent for systems with long-range interactions, 
we must formulate the problem in the microcanonical ensemble which is the fundamental one. 
We note that the terms involving the derivative of the density of states with respect to energy 
can be split in two parts according to 







9(E) dE 



9(E)Pj 



dP- 



(12) 



The terms with the E derivative would not have emerged if we had started from the Gibbs 
canonical distribution [36J. The equivalent hierarchy of equations for material particles in 
interaction is given in 
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2.2 Thermodynamic limit and mean field approximation 



Since systems with long-range interactions are generically spatially inhomogeneous, the usual 
thermodynamic limit N, V — > +00 with N/V fixed is clearly irrelevant. We define the proper 
thermodynamic limit of the point vortex gas as N — > +00 in such a way that the dimensionless 
temperature rj = f3N^ 2 and the dimensionless energy e = E/(N 2 ~f 2 ) are fixed. It is convenient 
to rescale the parameters such that 7 ~ l/N, E ~ 1, /3 ~ N and V ~ 1. Then, the total 
circulation T = Nj remains of order unity. We note that the ratio of dPj/dE on /3Pj is of 
order l/(Ej3) = l/(erjN). Therefore, in the thermodynamic limit N — > +00 with e, 77 fixed, the 
second term in Eq. (]12p is always negligible with respect to the first. Using this simplification 
in the second equation of the equilibrium hierarchy, we get 

^( ri ) = -(3(N - 1)7 2 / P 2 (r u r 2 )^dr 2 - (N - 1) 7 2 / ^^dr 2 , (13) 
or 1 J or 1 J or 1 oE 



OP 



(ri, r 2 ) = -/tfPfa, r 2 )^ - /3(N - 2) 7 2 / P 3 (r l5 r 2 , r 3 )^dr 3 . (14) 
or 1 or 1 J or 1 

We now decompose the two- and three-body distribution functions in the suggestive form 

P 2 (n, r 2 ) = Pi(n)Pi(r 2 ) + P 2 (rx, r 2 ), (15) 

P 3 (n, r 2 , r 3 ) = Pi(ri)Pi(r 2 )Pi(r 3 ) + ^(n, r 2 )Pi(r 3 ) 
+P 2 '(r l7 r 3 )P x (r 2 ) + P 2 '(r 2 , r 3 )P 1 (r 1 ) + i*(n, r 2 , r 3 ). (16) 

This decomposition is the counterpart of the first terms of the Mayer expansion in plasma 
physics. The P- are called the cumulants or the correlation functions. Inserting these de- 
compositions in Eqs. (IT3"l) - (IT4"j) . we find after simplification that the first two equations of the 
equilibrium hierarchy can be written 



dPi, , 



-/3(iV-l) 7 2 Pi(r 1 ) / P 1 {r 2 )^dr 2 
J or 1 



dri 

-W - 1)7 2 / %(T l7 r 2 )^dr 2 - (N - 1) 7 2 J ^^1, r 2 )dr 2 , (17) 



9ri 



(n,r a ) - (N-ltfPfa) f ^^( ri) r s )dr s = 

or 1 or 1 



+/3 7 2 A(ri)Pi(r 2 ) |^Hp 1 (r 3 )dr3-/5(iV-2)7 2 P 2 '(r 1 ,r 2 ) J Pl (r 3 )^dr 3 

WA(r 2 ) / ^P 2 '(r 1 ,r 3 )rfr 3 -/3(iV-2)7 2 Pi(r 1 ) / P^(r 2 ,r 3 )^r 3 
J or 1 J drt 

-P{N - 2) 7 2 J ^(n, r 2 , r 3 )^dr 3 , (18) 

where we have used Eq. (1171) to simplify some terms in Eq. (1TB]) . In the thermodynamic limit 
defined previously, it can be shown that the correlation functions P' n are of order jV-(»-i) [18]. 
Here, we shall just establish this result for the two-body distribution function P 2 assuming that 
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it holds at higher orders. We thus neglect the term P3, of order N~ 2 , in Eq. (|18p . This is the 
counterpart of the Kirkwood approximation in plasma physics. Then, considering the scaling 
of the terms in Eq. ([TBI) , we see that Pi ~ 1 and P 2 ' ~ (3^ 2 = rj/N = 0(1/N). Therefore, 

P 2 ( ri , r 2 ) = P 1 (r 1 )P 1 (r 2 ) + 0(1/N), (19) 

so that the mean field approximation p2(ri,r2) ~ Pi(ri)Pi(r2) is exact for iV — > +00. The 
coupling constant (3^ 2 ~ 1/N, scaling like the inverse of the point vortex number, plays a role 
similar to the "plasma parameter" in plasma physics. 



2.3 The mean field equilibrium distribution 

Taking the limit N —* +00 and using Eq. (1191) . the first equation (TTTj) of the equilibrium 
hierarchy becomes 

Vc^(r) = -/3 7 cj(r)V J w(r>(|r - r'|)dr', (20) 

where uj(t) = A r 7 P 1 (r) is the smooth vorticity field. After integration, this can be written in 
the form of the Boltzmann distribution 

w(r) = Ae-^ r \ (21) 

where 

V>(r) = J u(r')u(\r-r'\)dr', (22) 

is the stream function produced by the smooth distribution of point vortices. Therefore, the 
equilibrium density profile of the point vortices is determined by an integrodifferential equation. 
For the usual potential of interaction, satisfying Am = —5, we find that the equilibrium vorticity 
profile is determined by the Boltzmann-Poisson equation 

- Aip = Ae-^ T) . (23) 

These results can also be obtained by maximizing the Boltzmann entropy at fixed circulation 
and energy in order to obtain the most probable distribution of point vortices at statistical 
equilibrium [T7]. These results can be generalized so as to take into account the conservation 
of the angular momentum. In that case, the stream function ip in the Boltzmann distribution 
is replaced by the relative stream function ip' — ip + (l/2)f2x,r 2 [22] where is a Lagrange 
multiplier associated with the conservation of the angular momentum (the conservation of the 
linear impulse can be dealt with similarly |38j). 



3 Kinetic equation from a BBGKY-like hierarchy 
3.1 The BBGKY-like hierarchy 

We now address the out-of-equilibrium problem by using a methodology similar to the previous 
one. Our aim is to derive a kinetic equation for the evolution of the smooth vorticity profile 
cu(r, t) of the point vortex gas and to see whether or not it converges to the statistical equilibrium 
state (I2TI) . Integrating the Liouville equation ([3]) on Tj + i,...,rN, it is simple to construct a 
hierarchy of equations for the reduced distributions. It has the form 

f + E E v (^ i f + w-j) E / nj + i^ if-^rdv^ = 0. (24) 
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This is the counterpart of the BBGKY hierarchy in plasma physics. The first two equations of 
this hierarchy are 

dP x 



Of 



(N - 1) A J V(2 - l)P 2 (n, r 2 )dr 2 = 0, (25) 



BP BP B 

' V(2 - 1)^1 + ( N - 2)^- / V(3 - l)P 3 (r 1; r a , r 3 )dr 3 + (1 <- 2) = 0. (26) 



Bt dri dri 

For brevity, we have not written the variable t in the distribution functions. Inserting the 
decomposition ffT5l) in Eq. ff25l) . we first obtain 

9Pi 



(TV 



BP f BP 
1)^T y V(2 - l)Pi(r 2 )dr 2 + (iV - 1)— J V(2 - l)i*(n, r 2 )dr 2 = 0. (27) 



Next, substituting the decomposition (fl5|) and (TT6T) in Eq. (|26|) and using (|2"T|) to simplify some 
terms, we obtain 

f + V(2^)f + V(2^ 1 )P l(r2 )|i( ri ) 

9 



-Pi(r 2 )77- / V(3 - lJP^rOPi^drg 
oti y 

J V(3 -> l)P 2 '( ri ,r3)Pi(r 2 )rfr 3 
+ (iV - 2) A y V (3 - l)P 2 '(n, r 2 )P 1 (r 3 )rfr 3 
+(JV - 2)^- y V(3 -> l)P 2 '(r 2 , r 3 )P 1 (r 1 )rfr 3 

+(JV - 2)— / V(3 -> l)P 3 '(n, r 2 , r 3 )dr 3 + (1 <-> 2) = 0. (28) 

The equations (1271) and (128jl are exact for all N but the hierarchy is not closed. We shall now 
consider the thermodynamic limit defined in Sec. 12.21 Based on the scaling of the terms in 
each equation of the hierarchy, we argue that there exists solutions of the whole BBGKY-like 
hierarchy such that the correlation functions P- scale like l/iV- 7-1 at any time. This implicitly 
assumes that the initial condition has no correlation, or that the initial correlations respect 
this scaling (if there are strong initial correlations, like vortex pairs, the kinetic theory will 
be different from the one developed in the sequel). Recalling that P\ ~ 1, P 2 ~ 1/N and 
|V(z — > j)\ ~ 7 ~ 1/N, we obtain at order 1/N: 

BP f B f 

(N-l)-± V(2 -> l)Pi(r 2 )rfr 2 + N— / V(2 -> l)P 2 '(n, r 2 )dr 2 = 0, (29) 
oy x J dri J 



BP X 
~Bt 



BP' 
— - + 

at 



V(2 - 1) - J V(3 - l)P!(r 3 



)dr 3 



Pl(r2 W (ri ' 



BP' f B f 

+N J^ / V ( 3 ~> Wfo)^ + N— / V(3 - l)i*(r 2 , r 3 )P 1 (r 1 )rfr 3 + (1 <- 2) = 0. (30) 
Oti J <9ri J 

The three-body correlation function can be neglected. If we introduce the notations u = N^Pi 
(smooth vorticity field) and g = iV^P^ (two-body correlation function), these equations can be 
rewritten 

9 



^7 + (V ITT- + —V 2 -> 1 W2 — 
OT OTi 7 2 OTx 

+7T- / V ( 3 - ^fa' r 3' *)— dr 3 + (1 ^ 2) = 0. (32) 

For brevity, we have used the abbreviations oj\ = u(ri,t) and 0J2 = Ct>(r 2 ,t). We have also 
introduced the mean velocity in ri created by all the vortices 

(V) 1 = J V(2^1)^r 2 , (33) 

and the fluctuating velocity created by point vortex 2 on point vortex 1: 

V(2 - 1) = V(2 - 1) - 1<V)l (34) 

These equations (]3ip - (|32p are exact at order 0(1/ N). They form therefore the right basis to 
develop a kinetic theory. 

(i) Collisionless regime: If we consider the limit N —* +00 (for a fixed time t), noting that 
P' 2 = 0(1/ N) — > 0, we find that the smooth vorticity field ou(r,t) of the point vortex gas is 
solution of the 2D Euler equation 

^ + (V)Vcu = 0, (V) = -z x (35) 

where the stream function ip(r,t) is given by Eq. (1221) with u;(r,t) instead of oo(r). The 2D 
Euler equation is valid when the correlations between point vortices can be neglected, i.e. 
^2( r i) r 2)*) = Pi(Ti,t)P2(r2,t), which is the case for N — > +00. The Euler equation describes 
the collisionless evolution of the point vortex gas up to a time of order Nto (where to is the 
dynamical time) at least. In practice, N > 1000 so that the domain of validity of the 2D 
Euler equation is huge. The Euler equation is the counterpart of the Vlasov equation in plasma 
physics and stellar dynamics. It can undergo a process of violent relaxation towards a Quasi 
Stationary State [3TJ [321 [331 El] as discussed in Sees. 14.21 and 14.31 

(ii) Collisional regime: If we want to describe the collisional evolution of the point vortex 
gas, we need to consider finite iV effects. Equations (j3~TT) - (!32l describe the evolution of the 
system on a timescale of order Ntp. The equation for the evolution of the smooth vorticity 
field is of the form 

|f + ^<V)Vw = CWM, (36) 

where CV is a "collision" term analogous to the one arising in the Boltzmann equation. In 
the present context, there are not real collisions between point vortices. The term on the right 
hand side of Eq. f[36|) is due to the development of correlations between vortices as time goes 
on. It is related to the two-body correlation function g(ri,r 2 ,t) which is determined in terms 
of the vorticity by Eq. ( J32~l) . Our aim is to obtain an expression for the collision term Ctv[o;] 
at the order 1/N. The difficulty with Eq. ( 1321) for the two-body correlation function is that 
it is an integrodifferential equation. The second term is an advective term, the third term is 
the source of the correlation and the fourth term takes into account collective effects. In this 
paper, we shall neglect the contribution of the integral in Eq. (132I) . Then, we get the coupled 
system 
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dg_ 
dt 



d d 



9 



— — = 0. 38 

7 7 



The integral that we have neglected contains "collective effects" that are taken into account 
in the approach of Dubin & O'Neil [20] • However, their study is restricted to axisymmetric 
flows and makes a Markovian approximation. These assumptions are necessary to use Laplace- 
Fourier transforms in order to solve the integro-differential equation (1321) . Alternatively, if we 
neglect collective effects, we can obtain a general kinetic equation in a closed form (1521) that is 
valid for flows that are not necessarily axisymmetric and that can take into account memory 
effects. This equation has interest in its own right (despite its limitations) because its structure 
bears a lot of physical significance. Before deriving this general equation, we shall first consider 
the case of axisymmetric flows and obtain a simple explicit kinetic equation valid for such flows 
when collective effects are neglected. 

3.2 The Markovian axisymmetric equation 

For an axisymmetric flow, the vorticity field and the two-point correlation function can be 
written as u = ou(r, t) and g = g(rx, r 2 , Q\ — 9 2 , t), and the mean velocity as (V) = (V)g(r, t)ee. 
On the other hand, the projection of V(2 — ► 1) in the direction of ri is 

K„(2^1)= 7 I|H (39) 

where u 12 = n(r 1 ,r 2 ,6 l 1 — 6 2 ) is symmetric in t\ and r 2 (see Appendix |A|) . In that case, Eqs. 
([37D-(138]) become 

2 1 d f + °° , f 2n du . L , 7J 
-Qf = ~7 — g^r J r 2« r 2 / —9{ri,r 2 ,(p,t)d(p, 



o 



dq ^/ v i do du ( 1 d 1 d \ , .to . , 

£ + {a[r " l) ~ a(r " t)] i = -TAn^-^) ^ t] ' (41) 

where we have set = Q\—Q 2 and where Q(r, t) = (V)g(r, t)/r is the angular velocity of the mean 
flow. Taking the Fourier transform of Eq. (|4ip with respect to <fi and introducing the notations 
d = (l/ri)d/dri — (l/r 2 )d/dr 2 , io\ = uj(r 1 ,t), uj 2 = uj(r 2) t) and Af2 = Q(ri,t) — Q(r 2 ,t), we 
obtain 

^tp. _|_ imAQg m = ——mu m dujiuj 2 . (42) 
dt 7 

The Fourier transform of the potential of interaction u is discussed in Appendix |A] where 
explicit examples are considered. In terms of the Fourier transform of the correlation function, 
the kinetic equation (HU1) can be rewritten 



dt 



1 d f »> 

-27T7 2 — — / r 2 dr 2 y^mu m lmg m (ri,r 2 ,t). (43) 

r i dr i Jo ~ 



We shall assume that Im^ m (r 1; r 2 , t) relaxes on a timescale which is much smaller than the 
timescale on which w(r,t) changes. This is the equivalent of the Bogoliubov hypothesis in 
plasma physics. If we ignore memory effects, we can integrate the first order differential equation 
by considering the last term as a constant. This yields 

ft A 

9m(ri,r 2 ,t) = - dT-mu m e~ im ^ T d^ 2 , (44) 
Jo 7 
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where we have assumed that no correlation is present initially: g(t — 0) = 0. Then, we can 
replace Img m (ri, r 2 , t) in Eq. (143!) by its value obtained for t — * +00, which reads 

7T 

Im<7 m (ri, r 2 , +00) = mu m 5(m AQ)duiU2. (45) 

7 

Substituting this relation in Eq. ( l4"3l) . we obtain the kinetic equation 

1 = 2vr 2 7 -— / r 2 rfr 2X n, r 2 )S fii - ft 2 -^-^ - -^-^ , 46 
where we have defined 

X(ri, r 2 ) = ^ |m|u m (ri, r 2 ) 2 . (47) 

Til 

For the potential of interaction (11581) . this function is given by Eq. f II 6 1 j) and we recover the 
kinetic equation obtained in 



du\ 7 1 d 



dt 4 t\ dri J 



r 2 c?r 2 In 



, 1 duji 1 dojo , 
8{n x - tt 2 ) ~ -^ijr 1 ■ 48 



This equation, which ignores collective effects, is the vortex analogue of the Landau equation 
in plasma physics. We can show [231 [29] that it conserves T, E, L, that it satisfies an H- 
theorem (S > 0) and that the Boltzmann distribution (160]) is a particular steady state, but 
not the only one (see [22] for more discussion). Collective effects can be taken into account by 
keeping the contribution of the last integral in Eq. f)32p . For axisymmetric flows, these terms 
could be evaluated at the price of complicated calculations in the complex plane similar to 
those performed in plasma physics to derived the Lenard-Balescu equation from the BBGKY 
hierarchy. It would be interesting to make this derivation although it will not be attempted 
in the present paper. This would certainly lead to the kinetic equation derived by Dubin 
& O'Neil [2D] from a quasilinear theory of the Klimontovich equation. As we shall see, the 
consideration of collective effects is equivalent to replacing the bare potential of interaction by 
an "effective potential" . The resulting kinetic equation remains of the form of Eq. (146!) with a 
modified function Xp{ r ii r 2)- Therefore, as far as the general structure of the kinetic equations 
is concerned, our simple treatment is of interest. Furthermore, it can be generalized to non 
axisymmetric flows as considered in the next section. Finally, since the diffusion coefficient 
in Eq. (1461) does not diverge (contrary to the 3D Landau equation in plasma physics), the 
Lenard-Balescu treatment of collective effects is not necessary in our case for a first analysis. 



3.3 The general non Markovian kinetic equation 

The above kinetic equations assume that the flow is axisymmetric and rely on the assumption 
that the correlation function relaxes much more rapidly than the vorticity field. The Marko- 
vian approximation is expected to be a good approximation in the limit iV — > +00 that we 
consider since the vorticity profile changes only on a timescale of order Nt^ (where tp is the 
dynamical time) or even larger. However, for systems with long-range interactions, there can be 
situations where the decorrelation time of the fluctuations is very long so that the Markovian 
approximation may not be completely justified (this is the case for self-gravitating systems). 
For comparison, and for sake of generality, it can be of interest to derive non-markovian kinetic 
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equations for point vortices. For an arbitrary flow, Eq. (|38p for the correlation function can be 
written 



-(ri,t)-(r 2 ,t), 

7 7 



(49) 



where we have denoted the advective term by L (Liouvillian operator). Solving formally this 
equation with the Green function 



G(t : t') = exp j-jT £(r)dr 



we obtain 



g( ri ,r 2 ,t) = - I drG(t,t-r) [v(2 - 1)/- + V(l - 2)/- 



(50) 



-{Tx,t- T)-(r 2 ,t- T). 

7 7 



(51) 



The Green function constructed with the smooth velocity field (V) means that, in order to 
evaluate the time integral in Eq. (IBTj) . we must move the coordinates rj(£ — r) of the point 
vortices with the mean field flow (V)(r, t), adopting a Lagrangian point of view. Thus, in 
evaluating the integral, the coordinates rj must be viewed as rj(£ — r), where r^t — r) = 
Ti(t) — C ds (V)(rj(£ — s),t — s)ds. Substituting this result in Eq. (1371) . we obtain 



(9^! A/" — 1 ._ _. duj 



V u (2 



d 

] dr\ 



_d_ 

ar i Jo 

V u (l 



dr \ dr 2 y*(2-> l,t)G(t,t-r) 



a; 



u;(ri,t - r) — (r 2 ,t 

7 



(52) 



This returns the general kinetic equation obtained by Chavanis [23] with the projection operator 
formalism (note that we can replace V^{2 — > 1, t) by V^(2 — > 1, t) in the first term of the r.h.s. 
of the equation since the fluctuations vanish in average). It slightly differs from the equation 
obtained in [23] by a term (N — 1)/N in the l.h.s. This new derivation of the kinetic equation 
( 152"]) from a systematic expansion of the solutions of the BBGKY hierarchy in powers of 1/N 
is valuable because the formalism is much simpler than the projection operator formalism and 
clearly shows which terms have been neglected in the derivation. It also clearly shows that the 
kinetic equation is valid at order 1/N so that it describes the system on a timescale of 
order Nt D . In (SJ [2H], we had obtained this estimate a posteriori. 



3.4 Summary of the different kinetic equations 

Let us briefly summarize the different kinetic equations that appeared in our analysis. When 
collective effects are ignored, the kinetic equation describing the evolution of the system as a 
whole at order 1/N is 



~dt 



V (V)V " = ^/ <rfr / rfri ^ (1 



x<JV^(1^0)^- + V 1/ (0^1)^- 

v ' dr v y ' dr\ 



to(r, t 



Q)G(t,t-r) 
r). 



r - (r^t 

7 



(53) 
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If we make a Markov approximation and extend the time integral to infinity, we obtain 



x { v " (1 - + v " (0 - 11 ^ } w(r - ^ (r » (54) 

As we have indicated, the Markov approximation is justified for N — > +00 because the timescale 
Ntu on which uj changes is long compared to the timescale r corr for which the integrand in 
Eq. (154j) has significant support. We do not assume that the decorrelation time is extremely 
short so that, in the time integral, the vorticity and the vorticity gradient must be evaluated 
at r(t — r) and — r) where now rj(t — r) = r»(t) — Jq ds (V)(rj(t — s),t)ds. On the 
other hand, for an axisymmetric evolution, using the relation H153[) and rj(t — r) = rj(t) and 
^(t - r) = 0i(i) - n(n(t),t)r, Eq. ([53]) takes the form 



r dr d6 l I rrxdnV^l -»• 0, *) 
Jo Jo Jo 



9a; _ 1__9 
9t r 9r 



' 1 (9 1 (9 \ 

x y r (l_ ,t-r)( —)u(r,t-T)-{n,t-r). (55) 

1 r or riorij 7 



The integral on 9\ can be performed using Eq. (J98l) . and we get 

27T7-— I dr I r\dr\ "S^ m 2 u 2 m {r, r\) cos(mAf2r) 
r 9r y Jo 



(9a; 
9*~ 



19 19 



r dr ri dri 



ou(r,t — r)uj{ri,t — r). (56) 



This equation can also be obtained from the approach of Sec. 13.21 by keeping memory effects in 
Eq. (jUJ). If we make a Markovian approximation u (r 1; t — r) ~u(ri,t) and u (r,t — r) ~ o>(r, t), 
and extend the time integration to +00 in Eq. ( I55i) . we obtain 

duo 1 9 /"+°° r 2 ^ ^+00 

~dt = ~r~dr r J dT J Ml J rridri ^( 1 - °'*) 

x K (i _^ 0, t - r) f ~ - i A") W ( r| t )^( ri) t ). (57) 

The integral on r and 9\ can be performed using Eq. (1981) . see [231 EH] for details, and we get 
the kinetic equation (1461) . If we make the approximation a>(ri,t — r) ~ o;(ri,t) and u>(r,t — 
t) ~ o;(r, t) but keep the time integration from to t, we obtain the equation derived in 
[29] incorporating a function M(t) which regularizes the delta function occuring in Eq. (|46p . 
Finally, in Appendix [BJ we propose a simple heuristic kinetic equation that may be of interest. 



3.5 Discussion 

These kinetic equations possess a lot of interesting properties. Let us first consider the Marko- 
vian axisymmetric equation (1461) . The collisional evolution of point vortices is truly due to long 
range interactions because the current in r is caused by "distant collisions" with vortices lo- 
cated in ri 7^ r that can be far away. This is different from plasma physics and stellar dynamics 
where the collisions are assumed to be local [39]. Therefore, the current occurs only in velocity 
space and is due to "collisions" involving particles at the same location but having different 
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velocities V\ ^ v (recall that the position r in the point vortex system plays the same role as 
the velocity v in the plasma system). Therefore, in the case of stellar systems and plasmas, the 
collisional term is determined as if the system were spatially homogeneous. For these systems, 
long-range interactions manifest themselves only as mean field effects in the advective term 
(Vlasov) of the kinetic equation (see Eq. (44) of |40j). By contrast, the point vortex gas is the 
first system where collisions involve distant particles. The collisional evolution in r is due to a 
condition of resonance Q(ri,t) = Q(r,t) with point vortices in n 7^ r that have the same an- 
gular velocity. Clearly, this condition can be satisfied only when the profile of angular velocity 
is non-monotonic [201 [23] . Therefore, the evolution stops when the profile of angular velocity 
becomes monotonic (so that there is no resonance) even if the system has not reached statistical 
equilibrium. In that case, the system settles on a Quasi Stationary State (QSS) that is not the 
Boltzmann distribution (|2T|) predicted by statistical mechanics [29] . On the timescale Ntp on 
which the kinetic theory is valid, the collisions tend to create a monotonic profile of angular 
velocity. Since the entropy increases monotonically, the vorticity profile tends to approach the 
Boltzmann distribution but does not attain it in general because of the absence of resonances. 
The Boltzmann distribution may be reached on longer timescales, larger than Ntjj. To de- 
scribe this regime, we need to determine terms of order N~ 2 or smaller in the expansion of the 
solutions of the BBGKY hierarchy for N — > +00. This implies in particular the determination 
of the three-body correlation function, which is a formidable task. It is interesting to note that 
the markovian axisymmetric kinetic equation (I46p conserves all the integral constraints of the 
point vortex dynamics (circulation, energy, angular momentum) and satisfies an H-theorem for 
the Boltzmann entropy, so that the entropy is non- decreasing S > However, as we have 
indicated previously, this kinetic equation does not in general converge towards the Boltzmann 
distribution. This is because this kinetic equation admits an infinite number of stationary solu- 
tions among which the Boltzmann distribution is just a particular case (see [29] for a detailed 
discussion). This is at variance with the Landau and Lenard-Balescu equations which always 
converge towards the Boltzmann distribution [40]. In these equations, the collisional evolution 
of the system is also due to a condition of resonance k ■ vi = k ■ v (see Eq. (49) of [40]) but the 
Boltzmann distribution is the only steady state of these kinetic equations. As noted in [29] , the 
kinetic theory of point vortices is more closely related to the kinetic theory of one-dimensional 
systems with long-range interactions (like the HMF model) for which the collision term vanishes 
identically at order l /N [4711421140]. 

Let us now consider the more general kinetic equation (|52|) . We can prove by a direct 
calculation that this equation conserves the angular momentum and the linear impulse (see 
Appendix D of [23]). The conservation of the energy is more difficult to establish by a direct 
calculation but since Eq. (I52p is exact at order 0(1/N), the energy must be conserved (the 
integral constraints must be conserved at any order). Finally, we note that we cannot prove the 
if-theorem. It is only when additional approximations are made (markovian approximation) 
that the if-theorem is obtained (see Sec. 13.21 and |29j). To be more precise, let us compute 
the rate of change of the Boltzmann entropy for point vortices S = — J - In -dr with respect 
to the general kinetic equation (152]) . After straightforward manipulations, it can be put in the 
form 



If If* 
S = — - / drdri / dr 

2Y J uuJi Jo 



^(,.0)^ + ^(0.,)^ 



xG(t,t- r) 



(58) 



We note that, because of memory terms, the monotonic increase of the entropy is not granted. 
In the case of point vortices, the decorrelation time is much shorter than the relaxation time 
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(of order Ntu or larger) so that the markovian approximation is justified for N — > +00. In that 
case the entropy monotonically increases as shown explicitly for axisymmetric flows (see Sec. 
13.21 and [29]). However, as we have already indicated for the axisymmetric markovian equation 
(I46p . even if S > and E = T = 0, this does not imply that the kinetic equation will relax 
towards the Boltzmann distribution of statistical equilibrium [2JJ]. Indeed, the relaxation can 
stop before in the absence of resonances. The same remark applies to the more general equation 
(I52p . valid for non-axisymmetric flows, although this is more difficult to see. In order to make 
it clearer, one possibility would be to use the timescale separation between the dynamical 
time to and the collisional time t co u ~ Nto and derive an "orbit-averaged" kinetic equation 
in terms of appropriate variables similar to the angle-action variables used in other contexts. 
In that case, we would get a simpler kinetic equation, similar to the one derived in [13], and 
exhibiting an appropriate form of "resonances" between different orbits. This would generalize 
the condition of resonance Q(r) = Q(r') associated to Eq. (j4"6l) to the case of non-axisymmetric 
flows. Note that a phenomenological equation, valid for general flows, and exhibiting a form of 
"resonances" required to ensure the conservation of the energy is provided by Eq. (137) of [23] 
(see also Appendix IB"]) . 



3.6 The growth of correlations 

In Sec. 13.21 we have derived the equation satisfied by the two-body correlation function 
g(ri,r 2 ,<f),t) at order 1/N for an axisymmetric evolution. To derive the kinetic equation PS"]) , 
we only had to determine the imaginary part of the Fourier transform of g for t — > +00. In 
this section, we discuss the growth of the two-body correlation function in more detail. It is 
determined by the equation (}4"T]) . i.e.: 

7^7. + An ^I = ^I duJlUj2 - ( 59 ) 

ot ocp 7 ocp 

We shall assume that the initial vorticity profile is the Boltzmann distribution of statistical 
equilibrium 

uj = Ae -^ +1 ^ r2 \ (60) 

This distribution is a stationary solution of Eq. (1461) . Therefore, the r.h.s. of Eq. (I5T?|) is 
independent on time. Using Eq. (|6"UI) and the relations 

(V) e (r,t) = -^-(r,t) = n(r,t)r, (61) 

we find that 

duiu 2 = P^ui^AQ. (62) 

Substituting this result in Eq. flo^l) and introducing the function h = h((p, r 1; r 2 , t) through the 
relation g = ujiuj 2 h, we get 

M + A % = V (63) 
Taking the Fourier transform of the foregoing equation and integrating on time, we obtain 

h m (t) = -i/3 [ drmAQu m e- imAnr = /% m (e" imAm - 1). (64) 
Jo 
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Therefore 



9m(ri, r 2 , t) = (3u m {e- %m ^ u - 1)^ 2 . (65) 

We note that Im(g m ) has a limit (I45p for £ — > +00, while Re(g m ) has no limit. Taking the 
inverse Fourier transform of Eq. (1651) . we obtain 



g(<t>, ri, r 2 , t) = (3 [u(r 1 ,r 2 , <p ~ Afit) - u(ri, r 2 , 0)] uw 2 . (66) 

This equation describes the growth of two-body correlations in an axisymmetric flow assuming 
that the one-body distribution is given by the Boltzmann distribution. For the potential of 
interaction (1 1 5 8 [) written in the form 

U12 = —— \n(r\ + r\ - 2r 1 r 2 cos0), (67) 

47T 

the correlation function (|66|) is explicitly given by 



g{<f>,n,r 2 ,t) = --phi 



r\ — 2r\r 2 cos(0 — Af2i) 

uj\ijj 2 . (6? 



r\ + r 2 — 2r\r 2 cos 1 



We note that g has an oscillatory behavior, so that it has no limit for t — > +00. We recall that 
our approach is valid at order 1/N. Therefore, it corresponds to a linear regime extending on 
a timescale of order Nto- The two-body correlation function may reach the value it has at 
statistical equilibrium on a longer timescale but next order terms in 1/N in the developement 
must be taken into account. 



4 Kinetic equations from a quasilinear theory 

In this section, we show that the general kinetic equation (1521) describing the collisional evolution 
of the point vortex gas can also be derived from a quasilinear theory of the Klimontovich 
equation. We will compare the results with the quasilinear theory of the 2D Euler equation 
developed in [5U] to describe the process of violent relaxation in the collisionless regime. 

4.1 The slow collisional evolution of point vortices 

The exact vorticity profile of a gas of point vortices is a sum of Dirac functions given by Eq. 
([T]). It satisfies the equation 

^ + u d VuJ d = 0, (69) 

where = — z x is the exact velocity field created by uj^. Equation (169]) is the counterpart 
of the Klimontovich equation in plasma physics. It should not be confused with the 2D Euler 
equation (1551) or (ISTT) . the counterpart of the Vlasov equation, which has the same mathematical 
structure but which applies to the smooth vorticity field uo. The 2D Euler equation is valid 
during the collisionless regime (see next section) while the Klimontovich equation is exact and 
strictly contains the same information as the Hamiltonian equations (T5]). We now decompose 
the exact vorticity field in the form Ud = uj + 5uj where uj = (oud) is the smooth vorticity and Sou 
the fluctuation around it. Substituting this decomposition in Eq. (1691) and locally averaging 
over the fluctuations, we get 

^ + Luj = - (5uV5lu) , (70) 
ot 
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where L = u ■ V is an advection operator constructed with the smooth velocity field (note that 
u coincides with the field (V) introduced previously). Subtracting Eq. (I70p from Eq. (|69j) and 
neglecting non linear terms in the fluctuations, we obtain the following linearized equation for 
the evolution of the fluctuations 

95ou . _ . 

-— + L5u = -5uVu. (71) 
at 

Equations (1701) and (ITT]) form the basis of the quasilinear theory of the point vortex gas. These 
equations have been studied by Dubin & O'Neil [20] in the case where the smooth vorticity 
profile is axisymmetric. In that case, Eq. (J7T|) can be solved exactly with the aid of Laplace- 
Fourier transforms and a kinetic equation which takes into account collective effects can be 
obtained. This is the counterpart of the quasilinear theory in plasma physics that is used 
to derive the Lenard-Balescu equation from the Klimontovich equation [H]. In the present 
work, we shall proceed differently so as to treat the case of systems that are not necessarily 
axisymmetric and not necessarily Markovian. Our method avoids the use of Laplace-Fourier 
transforms and remains in physical space. This yields expressions with a clear interpretation 
which enlightens the basic physics. The drawback of our approach, however, is that it neglects 
collective effects. The formal solution of Eq. (!7T|) is 

5u(t) = G(t, 0)5oj(0) - f drG(t : t-r)5u(t-r)Voj(t-r), (72) 

Jo 

where G is the Green function associated with the advection operator L corresponding to the 
smooth mean field and we have noted u(t) = u)(r,t) and 5u(t) = 5u(r,t) for brevity. On the 
other hand, the perturbation of the velocity field is related to the perturbation of the vorticity 
through 

Suit) = - ( V(l -> 0)6u 1 (t)dr 1 , (73) 
7 J 

where refers to the position r. Therefore, considering Eqs. (1721) and (!73|) . we see that the 
velocity fluctuation 8u(t) is given by an iterative process: 5u(t) depends on 5ui(t) which itself 
depends on Sui(t — r) etc. We shall solve this problem perturbatively @ in the thermodynamic 



1 As shown in Sec. 12. 2\ the proper thermodynamic limit corresponds to N — > +oo in such a way that the 
individual circulation 7 ~ 1/N and the domain area V ~ 1. This implies that |r| ~ 1. We also have cu ~ 1 and 
5u> ~ 1/VN so that |u| ~ 1 and 5\u\ ~ l/\/N. With these scalings, we see that the terms that we have kept in 
Eq. ([71]) are of order u8u> ~ 1/y/N and u>Su ~ 1/y/N while the nonlinear terms that we have neglected are of 
order SujSu ~ 1/N <C 1/yN. We also note that the l.h.s. of Eq. ([70]) is of order u> ~ 1 while the r.h.s. of Eq. 
([70| is of order Sui6u ~ 1/N. Then Eq. ([70]) can be rewritten dtu> + Loj = (l/N)C(u) where the advective term 
is of order 0(1) and the collision term is of order 1/N. Therefore, this equation describes the evolution of the 
system on a timescale ~ Ntp. For N — * +00, it reduces to the 2D Eulcr equation dtuj + Loj = 0. In conclusion, 
the quasilinear theory developed in this section is equivalent to the truncation of the BBGKY hierarchy at the 
order 1/N. It amounts to neglecting three-body and higher correlations. 

2 It is at that stage of the dcvclopcment that we neglect some collective effects. In the approach of Dubin & 
O'Neil [20] for axisymmetric flows, Eqs. ([72]) and ([73]) can be solved exactly. 
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limit N — > +00 defined in Sec. 12. 21 To order 1/N we get 

(6uV6u) = f drtVCL 0)G 1 (t,0)G(t,0)(6u 1 (0)6uj(0)) 

7 or^ J 

1 d r* r 

dr / dTidT 2 V(l -> 0)Gi(t, * - r)G(t, * - r) 



x |y(2 - 0)(^(t - r)^ 2 (t - r))J^(t - r) 

+V"(2 - 1)<M* - r)5u 2 (t - r))||(t - r)}. (74) 
Now, the fluctuation is exactly defined by 

6u(t, t) = Yl - r *W) " w ( r > *)• ( 75 ) 

i 

Therefore, we obtain 

(Surfuz) = (^7 2 5( ri - r,)5(r 2 - r,-)) + (^7^^ - r,)^ - r< )) 

~(y^(ri - r*)^) - (^7^(r 2 - + ^1^2- (76) 

* 3 

To evaluate the correlation functions, we average with respect to the smooth distribution 
LOi/(N^) or UiUjj I \N^) 2 . This operation leads to 



N - 1 

(<5cji<k; 2 ) = — — — u^ 2 + 7^i<5(ri - r 2 ) - u) x uj 2 - uj 2 u 1 + cu^, (77) 



so that, finally, 



(Surfu^) = 7Wi<J(r 1 - r 2 ) - —u^jj 2 . (78) 



Substituting this result in Eq. (1741) . we obtain 

W = - 0)>^+ ^jfdr I drxV^l - 0)G(M~t) 

V(l - 0) Wl (t - r)^(t - r) + V(0 - - r)||(t - r)}, (79) 

where we have regrouped the two Greenians G and G\ in a single notation for brevity. Finally, 
replacing this expression in Eq. (]70p . we obtain the kinetic equation 

du N - 1 d f* f V M 

_ + _ uVw = 7 _ jf *■ y *,_(! o)G(t, < - r) 

^0)^,0^} . (80, 



X 



t-T 



This is identical to the general kinetic equation ()52[) obtained from the projection operator 
formalism or from the BBGKY-like hierarchy. We note that the term of order 1/N in the l.h.s. 
comes from the first term in Eq. (1721) . It corresponds to the mere advection of the vorticity 
fluctuation by the smooth velocity field in Eq. ([71]) . i.e. ignoring the coupling between the 
velocity fluctuations and the smooth vorticity (r.h.s. of Eq. (1711) ) which gives rise to the 
collision term. 
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4.2 The violent collisionless evolution of point vortices 

To leading order in N — > +00, the smooth vorticity profile of the point vortex gas is solution 
of the 2D Euler-Poisson system 

^ + uVcj = 0, cu = -A^. (81) 

This is the counterpart of the Vlasov-Poisson system in stellar dynamics and plasma physics. 
The 2D Euler equation describes the collisionless evolution of the point vortices due to mean 
field effects before collisions come into play on a timescale Nto or larger. Starting from an initial 
condition which is dynamically unstable, the 2D Euler-Poisson system develops an intricate 
filamentation at smaller and smaller scales. In this sense, the fine-grained vorticity u)(r,t) 
never achieves equilibrium. However, if we locally average over the filaments, the resulting 
"coarse-grained" vorticity uJ(r, t) will achieve a steady state on a timescale ~ to [SJ. Since the 
2D Euler equation is only valid in the collisionless regime t <C t co u, this corresponds to a quasi- 
stationary state (QSS) that will slowly evolve due to the effect of collisions occuring on a longer 
timescale ~ Ntp or larger. We can try to predict this QSS in terms of a statistical mechanics 
of the 2D Euler equation, using the approach of Miller- Robert-Sommeria [32l [33] IE This is the 
2D hydrodynamic version of the theory of violent relaxation proposed by Lynden-Bell [31] for 
collisionless stellar systems based on the Vlasov equation. In the case where the fine-grained 
vorticity w(r,t) takes only two values and ctq, the statistical equilibrium state maximizes the 
mixing entropy 

W = - ({?- in * + (l - ^ Infl - °) )dr, (82) 



at fixed circulation and energy. This leads to the coarse-grained vorticity 

u = -, vV^- ( 83 ) 

Note that the mixing entropy (182]) is formally similar to the Fermi-Dirac entropy and the 
equilibrium distribution (183"]) is formally similar to the Fermi-Dirac distribution. An effective 
"exclusion principle" , similar to the Pauli principle in quantum mechanics, arises in the theory 
of violent relaxation of continuous vorticity fields because the different levels of vorticity cannot 
overlap. We also note that the violent relaxation of point vortices does not lead to a segregation 
of the vortices according to their circulation (in the multi-species case) because the circulation 
of the individual vortices does not appear in the 2D Euler equation. Finally, we stress that 
the MRS theory is based on an assumption of ergodicity. Indeed, it implicitly assumes that 
the vorticity mixes well so that the QSS is the most mixed state compatible with the integral 
constraints of the 2D Euler equation. This may not always be the case as discussed in Sec. 14.31 
A kinetic theory of the process of violent relaxation has been developed in [3D] with the aim 
to determine the dynamical equation satisfied by the coarse-grained vorticity field lJ(t, t). This 
approach is based on a quasilinear approximation of the 2D Euler equation that is formally 
similar to that developed in the previous section (but with a completely different interpretation). 
In Sec. 14.11 the subdynamics was played by (a sum of ^-functions) and the macrodynamics 
by uj (a smooth field). The smooth field averages over the positions of the ^-functions (point 
vortices) that strongly fluctuate. In the phase of violent relaxation, the "smooth" field uo 



3 In these works, the 2D Euler equation is justified as a limit of the Navier-Stokes equation for inviscid fluids 
v — > 0. Since the 2D Euler equation also describes the collisionless regime of the point vortex gas, their approach 
can be applied in that context. 
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develops itself a finely striated structure and strongly fluctuates. Therefore, it is not smooth 
at a higher scale of resolution and a second smoothing procedure (coarse-graining) must be 
introduced. In that case, the subdynamics is played by u and the macrodynamics by To. The 
coarse-grained field averages over the positions of the filaments. The quasilinear theory leads 
to a kinetic equation for the coarse-grained vorticity of the form 



x 



§ + uW = e 2 ± jT dr J dr^Cl - 0)G(t, t - r) 
-(^Oh(, -^^ + ^(0-l)^o-^)^[ , (84) 



where e is the coarse-graining mesh size and u the velocity field produced by the coarse-grained 
vorticity (recall also that the ratio V(l — *■ 0)/j only depends on 17 and r according to Eq. 
flS}). This equation is expected to describe the late quiescent stages of the relaxation process 
when the fluctuations have weaken so that the quasilinear approximation can be implemented. 
It does not describe the early, very chaotic, process of violent relaxation driven by the strong 
fluctuations of the stream function. The quasilinear theory of the 2D Euler equation is therefore 
a theory of "quiescent" collisionless relaxation. This is the counterpart of the quasilinear theory 
of the Vlasov-Poisson system developed for collisionless stellar systems JUJUBE]. 

Equation (1841) is very similar, in structure, to Eq. (I80p for the collisional evolution of point 
vortices, with nevertheless three important differences: (i) the fluctuating velocity V(l — > 0) 
is replaced by the direct velocity V(l — > 0) because the fluctuations are taken into account 
differently, (ii) The vorticity u in the collisional term of equation flHDj) is replaced by the 
product EJ(o"o — uj) in equation ( 184]) . This nonlinear term arises from the effective "exclusion 
principle" accounting for the non-overlapping of vortex patches in the collisionless regime. This 
is consistent with the Fermi-Dirac-like entropy (1821) and Fermi-Dirac-like distribution (155]) at 
statistical equilibrium, (iii) Considering the dilute limit uj <C cr to fix the ideas, we see that the 
equations (184]) and (IBUj) have the same mathematical form differing only in the prefactors: the 
circulation 7 of a point vortex in Eq. (1801) is replaced by the circulation o"oe 2 of a completely 
filled macrocell in Eq. (I84p . This implies that the timescales of collisional and collisionless 
"relaxation" are in the ratio 

tncoll 7 / o r \ 

J. (85) 

tcoii o- e z 

Since o"oe 2 3> 7, this ratio is in general quite small implying that the collisionless relaxation 
is much more rapid than the collisional relaxation. Typically, t nco u is of the order of a few 
dynamical times to (its precise value depends on the size of the mesh) while t co u is of order 
Ntr>, or larger. The kinetic equation (l84"j) conserves the circulation, the angular momentum 
and, presumably, the energy. By contrast, we cannot prove an if-theorem for the MRS entropy 
( 182]) . Indeed, the time variation of the MRS entropy is of the form 

Smrs = Je 2 / drdr 1 - — ^ — [ drQ(t)G(t,t - r)Q(t - r), (86) 



Q(t) 



— (1 Ojwi^o - u x )— + —(0 lV(a - UJ)^ 



(87) 



and its sign is not necessarily positive. This depends on the importance of memory effects. 



21 



4.3 Discussion: incomplete violent relaxation 



Even if Eq. (1841) conserves the energy and the circulation and increases the MRS entropy (1521 
monotonically, this does not necessarily imply that the system will converge towards the MRS 
distribution There are many cases where the MRS theory provides a good prediction of 

the QSS [4"5l |4"§] . However, it has also been observed in some experiments [5U] and numerical 
simulations [51] that the QSS does not exactly coincide with the strict statistical equilibrium 
state predicted by the MRS theory because of the complicated problem of incomplete relaxation 
[52J. This is usually explained by a lack of ergodicity or "incomplete mixing". Here, we try to 
be a little more precise by using the kinetic theory. There can be several reasons of incomplete 
relaxation: 

(i) Absence of resonances: Very few is known concerning kinetic equations of the form of 
Eq. (1841) and it is not clear whether the MRS distribution (!83|) is a stationary solution of that 
equation (and whether it is the only one). In order to simplify Eq. (I84|) . we shall assume 
that the coarse-grained vorticity field is axisymmetric and that the correlations relax on a 
timescale which is much shorter than the typical time on which the coarse-grained vorticity field 
changes (Markovian approximation). Although this approximation was justified to describe 
the collisional relaxation of point vortices (because of the timescale separation between the 
dynamical time and the collision time), this approximation is not clear for the process of 
violent collisionless relaxation where memory terms can be important. However, with this 
approximation we can re-do the calculations of the previous sections and obtain a kinetic 
equation of the form 



duJ 
~dt 



4r dr 



+00 



r±dri In 



r> 



Qi) <{ -ui(o- — uji) 



dui 
dr 



1 dui 
—uj(a -u)—— 
r\ ori 



This equation conserves the circulation, the energy and the angular momentum and satisfies 
an H-theorem for the mixing entropy (I82"j) . However, as already discussed for the collisional 
evolution of the point vortices, it does not relax in general towards the statistical equilibrium 
state, here the MRS distribution fl83|) because of the absence of resonances. The same conclusion 
probably applies to the more general equation (|84"|) and to the heuristic equation (148) of [23j 
which also conserve E, T, L and increase the mixing entropy ( 182|) but do not obligatory relax 
towards the MRS distribution ( l83i) . The system tries to approach the statistical equilibrium 
state (as indicated by the increase of the entropy) but may be trapped in a QSS that is different 
from the statistical prediction (!83|) . This QSS is a steady solution of Eq. (1841) which cancels 
individually the advective term (l.h.s.) and the effective collision term (r.h.s.). This determines 
a subclass of steady states of the 2D Euler equation (cancellation of the l.h.s.) such that the 
complicated "turbulent" current J in the r.h.s. vanishes. This offers a large class of possible 
steady state solutions that can explain the deviation between the QSS and the MRS statistical 
equilibrium state fl83l) observed, in certain cases, in simulations and experiments of violent 
relaxation. One may argue that nonlinear terms are needed in the kinetic theory in order to 
obtain an equation that relaxes towards the statistical equilibrium distribution fl83l) . In that 
case, we must develop a kinetic theory that goes beyond the quasilinear approximation. 

(ii) Incomplete relaxation in space: The turbulent current J in Eq. (1541) is driven by the 
fluctuations u 2 = Cj 2 = u 2 — ZJ 2 of the vorticity (generating the fluctuations 5u of the velocity) 
|30j . In the "mixing region" where the fluctuations are strong, the vorticity tends to reach 
the MRS distribution fl83|) . As we depart from the "mixing region", the fluctuations decay 
(co>2 — ► 0) and the mixing is less and less efficient. In these regions, the system takes a long 
time to reach the MRS distribution (1831) and, in practice, cannot attain it in the time available 
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(see (iii)). In the two levels case, we have u>2 = TD(ctq ~~ w). Therefore, the regions where uj — ► 
or uJ — > (To do not mix well (the diffusion current J is weak) and the observed vorticity can be 
sensibly different from the MRS distribution in these regions. This concerns essentially the core 
(ZU — > Co) and the tail (uJ — > 0) of the vorticity distribution. This result, derived from the kinetic 
theory, is consistent with what is observed in experiments [5UIET]: the core vorticity decreases 
less than predicted by the MRS statistical theory while the tail of the vorticity profile decreases 
more rapidly. For the same reason, it can also explain why the vorticity peak is remarkably well 
conserved during a merging process as observed in 2D decaying turbulence [53J. 

(iii) Incomplete relaxation in time: during violent relaxation, the system tends to approach 
the statistical equilibrium state ( 1531) . However, as it approaches equilibrium, the fluctuations 
of the velocity field, which are the engine of the evolution, become less and less effective to 
drive the relaxation. This is because the scale of the fluctuations becomes smaller and smaller 
as time goes on. This effect can be taken into account in the kinetic theory by considering 
that the correlation length e(t) decreases with time so that, in the kinetic equation (151)) . the 
prefactor e(t) — > rapidly. As a result, the "turbulent" current J in Eq. (I84p can vanish before 
the system has reached the statistical equilibrium state (1531 . In that case, the system can be 
trapped in a QSS that is a steady solution of the 2D Euler equation different from the statistical 
prediction (1531) . 

5 Relaxation of a test vortex in a bath 

In this section, we study the relaxation of a test vortex in a bath of field vortices. Specifically, we 
consider a collection of N point vortices at statistical equilibrium (thermal bath) and introduce 
a new vortex in the system. To leading order in N — > +oo, the point vortex is advected by 
the mean flow. However, due to finite N effects, the test vortex undergoes discrete interactions 
with the vortices of the bath and progressively acquires their distribution. We wish to study 
this stochastic process. The probability density P(r,t) of finding the test vortex in r at time t 
is governed by a Fokker-Planck equation involving a term of diffusion and a term of drift. In 
our previous papers, we obtained these terms from a linear response theory [51] or from the 
projection operator formalism (23]- In the present work, we obtain these terms directly from the 
equations of motion and show how collective effects can be included in the theory in the case 
of axisymmetric flows. Our approach is also valid if the bath is made of a distribution of field 
vortices that evolves slowly, so that it can be assumed stationary on a timescale (N/hx N)to 
which is the typical relaxation time of the test vortex in the bath (see [29] for details). 

5.1 Diffusion coefficient 

The increment of the position of the test vortex between t and t — s due to the fluctuations of 
the velocity is 



After standard calculations (see, e.g., Sec. 4.2 of [SSj), the second moment of the increments 
of position can be rewritten in the form 




(89) 




(90) 
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We shall assume that the correlation function decreases more rapidly than r 1 . Then, taking 
the limit s — * +00, we find that the diffusion coefficient is given by the Kubo formula 



-00 



= ( = lim ( ) = / (V»(f)V(t - r))dr. 



2At / \ 2s 

On the other hand, after straightforward calculations (see, e.g., Sec. 4.1 of [55]), we obtain 
(V^V^t - r)) = JV(V*(1 -> 0, *)V(1 -> 0, t - r)) 
= [ dnV^l -»0,t)V(l ->0,t-r)-(r 

7 7 

Therefore, combining Eqs. (19T1) and fl92|) . we get 



(92) 



C+oo 

D^= drdr 1 V ,i (1^0,t)V i/ (1^0,t-r)-(r 1 ). (93) 
Jo 7 



For an axisymmetric system, the diffusion coefficient is given by 

= *„ /(^ 

2At / ^+oo\ 2s 







leading to 

^+00 p2n /*+oo 

D= / dr d6t ridriV r (l->0,t)V r (l->0,t-r)-(ri). (95) 
Jo jo Jo 7 

If we neglect collective effects, the velocity created by a field vortex 1 on the test vortex is 
given by (see Appendix 1X1) : 

V r (l -»• 0, t) = 17- ^ mv iro ( Wl ». (96) 

m 

At time t — t, we have 

K(l -»• 0,t -r) = i7-^m« m e !m(9(< " r) " fll(t " T)) . (97) 



r 

m 



To leading order in A" — > +00, the point vortices are advected by the mean field velocity so 
that 9i(t — r) = 6i — fi( r «) r an d r i(^ — T ) = r i where = Tj(t) and 0j = 6>j(t) denote their 
position at time t. Thus, we get 



V r {l -> 0,t - r) = i 7 ~ ^mti ra e ,mW - A!lT) . (9* 



r 



Substituting this expression in Eq. (193]) and carrying out the integrations on 6 1 ! and r, we 
obtain after straightforward calculations 

2tt 2 7 f +OD 

D = ^^ ndrix(r,n)5(A^)a;(ri), (99) 
r Jo 

where the function x(r, ri) is defined in Eq. (1471) . If the profile of angular velocity is monotonic, 
we can use 5(AQ) = 8(r — ri)/\Q'(r)\ and we find that 

D(r) = 2At^|^(0, (100) 
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where S = rQ'(r) is the local shear. For the potential f)158p . we have 

1 + °° 1 1 

x(r,r) = —J2- = ^^^ (101) 

where In A = X^m=i m * s a logarithmically diverging Coulomb factor that has to be regularized. 
This leads to the following expression of the diffusion coefficient 

D{r) = 4 lnA M^ (r) - (102) 

This expression, with the shear reduction, was derived in Chavanis [SU [23] and Dubin & Jin 
[56] from the Kubo formula. Because of the divergence of the Coulomb factor, the value of the 
diffusion coefficient is dominated by the contribution of field vortices at radial distance r\ = r, 
justifying the local approximation made in [51]. Therefore, Eq. (11021) gives the dominant term 
in the diffusion coefficient ( J99l) of a test vortex even if the profile of angular velocity of the field 
vortices is non-monotonic. In practice, the Coulomb factor has to be regularized as discussed 
in detail in [56], [29] . It is then found that In A scales with the number of particles like \ In N in 
the thermodynamic limit N — > +oo, in agreement with the rough estimates in [54] . 



5.2 Drift coefficient 

In addition to its diffusive motion, a test vortex immersed in a bath of field vortices with spa- 
tially inhomogeneous vorticity distribution undergoes a systematic drift. The drift corresponds 
to the response of the field vortices to the perturbation caused by the test vortex, as in a 
polarization process. The test vortex modifies the density distribution of the field vortices and 
the retroaction of this perturbation on the test vortex causes its drift. The expression of the 
drift can be derived from a linear response theory starting from the Liouville equation as done 
in [53]. In this section, we show that it can also be obtained from the Klimontovich equation. 
This will make a close connection with the quasilinear theory developed in Sec. [31 

The introduction of a test vortex in a bath of field vortices modifies the vorticity profile 
cu(r, t) of the bath. Since this perturbation is small, it can be described by the linearized Euler 
equation 

^ + L5uo = -5uVu, (103) 
at 

whose formal solution is 

6u(t) = - [ drG{t,t-T)5u{t-T)Vu(t-T). (104) 
Jo 

We have used the fact that, initially, d~u(0) = 0. On the other hand, the perturbation of the 
velocity field in r is given by 

8u{r, t) = - J V(l -> 0)5u 1 (t)dr 1 + J V(l -> 0)5( ri - r P (t))dn. (105) 

The second term is the velocity created by the test vortex at position rp(t) and the first term 
is the fluctuation of the velocity due to the perturbed density distribution of the field vortices. 
Substituting Eq. (TTMj) in Eq. ffTUoT) . we obtain 

8u(r, t) = -- [ dr [ dnVfl -> 0)GAt, t - T)5u"At - r)^(t - r) 
7 Jo J or" 

+ f V(l -> 0)5(ri - r P (t))dn. (106) 
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This is an integral equation for 5u(r,t). For an axisymmetric flow, this equation can be solved 
exactly by using Laplace- Fourier transforms as done in Schecter & Dubin [5?]. In order to 
treat more general flows, we shall make an approximation which amounts to neglecting some 
collective effects. We solve Eq. (I106P by an iterative process: we first neglect the first term in 
Eq. (jl06p keeping only the contribution of the test particle. Then, we substitute this value in 
the first term of the r.h.s of Eq. (11061) . This operation gives 

8u(r, *) = ™ f Q dT J dridr 2 V(l _> 0)Gi(i, t - t)V v {2 -> 1) 
d ±(t - r)5(r 2 - r P (t ~r)) + J V(l - 0)£(n - r P (t))dn. (107) 



x 



This quantity represents the fluctuation of the velocity field in r caused by the introduction 
of a test vortex in the system and taking into account the retroaction of the field vortices. 
If we evaluate this expression at the position rp of the test vortex and subtract the second 
term (self-interaction), we obtain the drift experienced by the test vortex in response to the 
perturbation that it caused. Denoting finally by the position of the test vortex, we find that 
its drift is given by 

Koi = ~- fdr [ dv 1 V»(1^0,t)V»(0^1,t-T)^(v 1 (t-T)). (108) 



7 Jo J dr\ 



For a thermal bath, where the distribution of the field vortices is given by u}{v-\) = Ae ^( Tl \ 
we obtain 



V £* = Pf dT f ^i^(1^0,t)V(0^1,t-r)-V^(r 1 (t-r)V(r 1 ), 



(109) 



where we have used u(ri(t — r)) = u(r\(t)) since u is a stationary solution of the 2D Euler 
equation. This is equivalent to the result of the linear response theory based on the Liouville 
equation [51] but it is obtained here in a simpler manner. We can also obtain this result in a 
slightly different way. We approximate Su(r, t) in Eq. fl 1 3 [) by the velocity V(P — > 0) created 
by the test vortex so that 

— - + L5u = -V(P -> 0)Vcu. (110) 
at 

This equation can be solved with Green functions yielding 

Su{t) = - [ drG(t,t-T)V(P ->0,t-T)Vu>(t-r). (Ill) 
Jo 

This represents the perturbation of the distribution of field vortices caused by the introduction 
of the test vortex in the system. This perturbation produces in turn a velocity which causes 
the drift of the test vortex (by retroaction). If we substitute Eq. (11 111) in the first part of Eq. 
( 11051) and evaluate this quantity at the position of the test vortex, we recover Eq. fll08j) for the 
drift (see also Appendix [C]) . 

If we now consider an axisymmetric distribution of field vortices, the expression of the drift 
becomes 

i rt r2iT r+co j 

v r p° l = — dr de x rit z ri y r (i_o,t)y ri (o-^i,t-r)- r (n). (112) 

7 Jo Jo Jo ar 
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Using the identity f 1 1 5 3 j) and taking the limit t — > +oo, we get 

i r+°° r 2n r+°° r 
V r pol = - dr d6 1 r 1 dr 1 -K(1^0,t)K(1^0,t-r)— (n). (113) 
1 Jo Jo Jo n dr 

This is a sort of generalized Kubo relation involving the gradient of the density profile instead 
of the density profile itself. The nice similarity in the expressions of the diffusion coefficient 
K95\) and drift term U13\) is worth mentioning. The integrals on 9\ and r can be evaluated in 
the same manner as in Sec. 15.11 and we obtain 

y P oi = 2^7 f + °° driX ( r)ri )j( An) ^( ri ). (114) 
r J Q dr 

Now, the drift of the test vortex is due not only to the polarization process but also to the 
variation of the diffusion coefficient with r. As a result, the complete expression of the drift is 

V dr lf t^/^:\ = dD + vrK (115) 



t At I dr 
From Eqs. (1991) and (I114p . we obtain 

V r = 2- 2 7 1 rndrMn) (~ r Q- r ~ -^rj X^r^ACl)-, (116) 

where we have used an integration by parts in Eq. ( 11 14ft . Expressions ( |99l and (11161) for 
the diffusion coefficient and the drift term can also be obtained directly from the Hamiltonian 
equations, by making a systematic expansion of the trajectory of the point vortices in powers 
of 1/N in the limit iV — > +oo as shown in Appendix C of |29j . 

For a thermal bath, corresponding to the case where the field vortices are at statistical 
equilibrium, the vorticity profile is given by the Boltzmann distribution 

u(n) = Ae-^' {ri \ (117) 

where ip' = ip + (l/2)f2^r 2 is the relative stream function. Then, we have 

= -p^M = pj^in) - n L )n. (us) 

dri dr± 

Substituting this relation in Eq. (11141) . using the (^-function to replace fi(ri) by fi(r), using 
dtfj' /dr = (— Q(r) + VLL)r and comparing the resulting expression with Eq. ( 1991) . we finally find 
that 

V** = -/3-yD^-. (119) 
dr 

The drift is perpendicular to the relative mean field velocity and the drift coefficient (mobility) 
is given by a sort of Einstein relation £ = Dfi^y. We note that the drift coefficient and the 
diffusion coefficient depend on the position and we recall that the temperature is negative in 
cases of physical interest |16j . We also emphasize that the Einstein relation is valid for the drift 
Vf° l due to the polarization only, not for the total drift ( 11161) . We do not have this subtlety for 
the usual Brownian motion where the diffusion coefficient is constant. 

If we now consider a bath with a monotonic profile of angular velocity, using the same 
arguments as in Sec. 15.11 we find that Eq. (11141) reduces to 
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For the potential (|158p . using Eq. (11011) . we find that 



Vf o1 = - In A 



1 duj 
|£(r)| dr 



r). 



'121) 



Due to the diverging factor In A ~ ^ In N, this expression also gives the dominant term of the 
drift in the case where the vorticity profile is non monotonic. Comparing with Eq. fU02p . we 
find that the drift velocity is related to the diffusion coefficient by the relation 



ypol = D 



dhiuj 
dr 



;i22) 



This expression generalizes Eq. (11 1 Q p for a bath that is out-of-equilibrium. 



5.3 Collective effects 

As explained previously, one specificity of our approach is to develop a formalism that allows to 
describe flows that are not necessarily axisymmetric. However, its main drawback is to ignore 
collective effects. In the case of axisymmetric flows, these collective effects can be taken into 
account as in the study of Schecter & Dubin [57]. In this section, we briefly indicate how the 
preceding results can be generalized to account for these collective effects. 

For axisymmetric flows, Eqs. (11031) and H 1 5 [) can be written (we consider here the usual 
situation where the potential of interaction between vortices is solution of the Poisson equation): 



35uj 



fl(r,t) 



dScu 1 88 ip 8u> 
86 r 89 dr 



0. 



(123) 



0* 

dr' 



+ 



18 Id 

+ 



r dr r 2 86 2 



5tjj = -8u - 7-<5(r - r P {t))8{9 - 9 P {t)), 



(124) 



where (rp, 6p) are the coordinates of the test vortex. These equations can be solved by taking 
the Laplace-Fourier transform of 5u and 5ip. Returning to physical space, the perturbed stream 
function can finally be written [57] : 



<ty(r, M)=7E e im(e - 6p) U m (r, r P ), 



'1251 



where 



U m (r, r P ) 



AtiH 



a+ic °G(r,r P ,m,Sj^ 



e ds, 



(126) 



where G is the Green function solution of 

,2 



d 2 1 d 



m 



im 



Iduj 



s + imVL{r) r dr 



G(r, rp, m, s) 



5{r — r P ) 



(127) 



Therefore, when we take into account collective effects, the radial velocity created by point 
vortex 1 on point vortex (say) is given by 



V r (l -> 0) = ^ S ^ j ime im{e - t>l) U m (r,r 1 ). 



128) 
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Neglecting collective effects amounts to neglecting the last term in brackets in Eq. (11271) . It 
then reduces to the usual Poisson equation where the vorticity field is due to a single point 
vortex. Then, G is equal to Gbare = — 27rw m (r, rp) so that U m {r, r±) is replaced by u m (r, ri) in 
Eq. (11281) . This returns the bare velocity (196]) created by point vortex 1 on point vortex 0. 

In the computation of the diffusion coefficient, we can take into account collective effects by 
replacing u m by U m in Eq. (19~B1) . This yields Eqs. (I9"9l and (1 1001) where x( r 5 r i) is replaced by 

X (r,r 1 ) = J2\ m \\Um(r,r 1 )\ 2 . (129) 

m 

In fact, for r — r%, the series diverges for large m indicating that the main contribution to 
the diffusion coefficient is due to close interactions, justifying a local approximation. We thus 
qualitativelty understand that collective effects will be negligible. For large m, we can replace 
\U m (r,r)\ 2 by \u m (r, r)| 2 returning the result (11021) . 

In the computation of the drift, we can take into account collective effects as follows. The 
velocity created in by the introduction of the test vortex is 



V r (P -> 0) = -J2 im U m (r,r P )e im{e - ep) . (130) 



The bare velocity due to the test vortex is 

V r (P -> 0) = - im u m (r, r P )e im{e ~ Bp) . (131) 

m 

If we subtract Eq. (1 1 3 1 f) from Eq. (I130p . we get the velocity created in by the perturbation of 
the distribution of the field vortices caused by the introduction of the test vortex. Evaluating 
this velocity at the position of the test vortex, we obtain the drift experienced by the test vortex 
due to the polarization process 

V r po1 = ^^2im{U m (r,r) - u m (r,r)). (132) 



This can also be written 



V r po1 = -- J2 m Im \u m {r, r)j . (133) 



The series diverges for m — > +oo. If we replace lm[U m (r,r)] by its asymptotic behaviour for 
large m [57], this returns Eq. 



Finally, we conclude that, concerning the evaluation of the diffusion coefficient and drift 
term, collective effects play a negligible role since these quantities are dominated by close 
interactions. This gives further justification to the approaches developed in Sees. 15.11 and 15.21 



5.4 The Fokker-Planck equation 

Assuming that the evolution is axisymmetric, the probability density P(r, t) = P(r, t) of finding 
the test vortex in r at time t is governed by a Fokker-Planck equation of the form 

dP 1 d 
dt 2r dr 

This Fokker-Planck approach assumes that the stochastic process is markovian which is a good 
approximation in our case, as we have already indicated. It also assumes that the higher order 



r 9_ ((Ml p 
dr 1 At 



ld_ 

r dr 



rP 



(Ar) 



(134) 
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moments of the increment of radial position Ar play a negligible role. This is indeed the case 
in the N — > +00 limit that we consider since they are of order 0(N~ 2 ) or smaller. At order 
0(N~ l ), we have found that the second (diffusion) and first (drift) moments of the radial 
increments of position of the test vortex are given by 

2Ai At dr h v ' 

with 

d= 2ttW rrfrafanMAtyuin), (136) 
r Jo 



= yp° l 



2vr 2 7 



00 dui 
drix(r,ri)5(Afi) — (n). 

dr 



(137) 



The Fokker-Planck equation (11341) can be written in the alternative form 



dP 
~dt 



ld_ 

r dr 



dP 

r I D — Prj 

or 



(138) 



The two expressions ( 11341) and (11381) have their own interest. The expression f ll34j) where 
the diffusion coefficient is placed after the second derivative d 2 (DP) involves the total drift 
y drift _ (Ar)/Ai and the expression H138[) where the diffusion coefficient is placed between 
the derivatives dDdP isolates the part of the drift rj = yv° lar due to the polarization. This 
alternative form f)138p has therefore a clear physical interpretation. Inserting the expressions 
and (I137p of the diffusion coefficient and drift term in Eq. (jl38j) . we obtain 



BP Id r +co (\ 8 

— = 2ttV- / ri dr lX (r, n)5(n - n x ) 
at r or J \r or 



]__d_ 

r\ dr\ 



P{r,t)u( ri ). 



(139) 



For a thermal bath, using Eq. (11191) . the Fokker-Planck equation ( 11381) can be written 



dP 
~dt 



ld_ 

r dr 



rD(r) 



dP_ + p <W\ 
dr dr J 



(140) 



where D(r) is given by Eq. (11361) . For a steady bath with a monotonic vorticity profile, using 
Eq. ( 11221) . the Fokker-Planck equation ( 11381) can be written 



dP 
~dt 



ld_ 

r dr 



rD(r) 



dP 

dr 



P- 



dhxou 



dr 



(141) 



where D(r) is given by Eq. (I102p . These Fokker-Planck equations have been studied in detail 
by Chavanis & Lemou [29J for different types of bath distribution. The distribution of the test 
vortex relaxes to the distribution of the bath on a typical timescale (N/ In N)to but the relax- 
ation process is very peculiar and differs from the usual exponential relaxation. In particular, 
the evolution of the front profile in the tail of the distribution is very slow (logarithmic) and 
the temporal correlation function (r(O)r(t)) decreases like hit/t (for a thermal bath). This is 
due to the rapid decay of the diffusion coefficient D(r), like in the HMF model [4T| . 

In our previous papers, we have obtained Eq. (I139[) directly from the projection operator 
formalism (see Sec. 4.1 of [29J). This amounts, in the kinetic equation (f4"6l) . to replacing u(r, t) 
by the distribution P(r,t) of the test vortex and u)(ri,t) by the static distribution u(ri) of the 
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field vortices. This procedure transforms an integrodifferential equation (]46p into a differential 
equation (I139p . Then, the expressions (I135p - (I137I) of the diffusion and drift terms were obtained 
by identifying Eq. (1 1 3 9 [) with the Fokker-Planck equation (I134p . In the present paper, we have 
proceeded the other way round by first determining the moments (I135p - (I137I) in Sees. 15.11 
and 15.21 then inserting them in the Fokker-Planck equation (11341) . Although this procedure 
may appear more logical in some sense, the other approach based on the projection operator 
formalism is more powerful because it allows one to obtain more general equations that are non 
markovian and that relax the hypothesis of axisymmetry as discussed in the next section. 



5.5 More general kinetic equations 

It is instructive to compare the Fokker-Planck equation (11391) with the more general equation 
obtained from the projection operator formalism 

+ (V) VP = -^J/tJ dnF(l - 0)G(t, t - r) 

x { V(l - 0) A + V (0 - 1) A}p( r , t - r)^{v x ). (142) 

This equation can be obtained from Eq. (153"P . by replacing u(r, t) by P(r,t) and u(ri,t) by 
u)(r-\_). This is a sort of generalized "Fokker-Planck" equation involving a term of "diffusion" 
and a term of "friction" . However, strictly speaking, Eq. (11421) is not a Fokker-Planck equation 
because it is non-Markovian. We also note that the "diffusion" term appears as a complicated 
time integral of the velocity correlation function involving P(r,t — r). This can be seen as 
a generalization of the Kubo formula fl9B"p . Similarly the "drift" is a generalization of the 
expression obtained in (11081) with a more complicated time integral. If we consider a thermal 
bath where the distribution of the field vortices is the Boltzmann distribution, we get 



^ + (V) VP = A jf dr J dv.V^l - 0)G(t, t - r) 



x |V(1 -> 0) ■ V - /3 7 V(0 -> 1) ■ W( ri ) |P(r, t - r)-(ri). (143) 

If we come back to Eq. (11421) . make a Markovian approximation and extend the time integral 
to infinity, we obtain 



dP d f +oa f 

~dt + (V) VP = d^J dT J driV " (1 ~* 1 

x I V(l - 0) A + V^(0 - 1) A\ P(P , tf-iv^ (144) 



T 



dr v dr" J 7 

where we recall that the coordinates appearing after the Greenian must be viewed as explicit 
functions of time, i.e. r,j = r^t — r). On the other hand, for an axisymmetric system, using 
the relation f[T53|) and n(t - r) = r^t) and 9$ - r) = 0j(t) - fi(ri(t))r, Eq. takes the 

simplest form 

dP Id f l f 27T f + °° 

~dt = ~r~dr r J dT I Ml J rr ^ V r( l ^ M) 

xV,(l - 0, t - r) f-|- - 1 A ")P(r, t - r)-(n). (145) 
\r or r\ ar\ J 7 
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If we make a Markovian approximation P(r,t — t) ~ P{r,t) and extend the time integral to 
infinity, we get 



dP 1 d r+°° r 27T r+oo 

~ot = rfr r J dT J Ml J rridri ^( 1 °'*) 

xVr{1 ^ 0) t-T)f~ - -^-)p(r,t)-( ri ). (146) 
\r or n or i J 7 

This is a Fokker-Planck equation which can be put in the form (11381) with a diffusion coefficient 

/*+oo /*27r r-\-oc 

D= dr dO l I ri d ri y r (l_>o,t)V r (l->0,t-T)-(ri), (147) 
7o Jo Jo 7 

and a drift term due to the polarization 

1 r+oa r-2-K r+oc i 

r/ = -/ dr / d0i / rdnK(1^0)V;.(1^0,t-r) — (n). (148) 
7 Jo Jo Jo « r 

These expressions agree with Eqs. ( 1951) and f II 1 3 j) obtained directly from the equations of 
motion. After integration on r and 9\, we recover the Fokker-Planck equation (11391) with the 
expressions ( 11361) and H13TI) of the diffusion coefficient and drift term. 



6 Conclusion 

In this paper, we have developed the kinetic theory of point vortices in two-dimensional hy- 
drodynamics initiated in [23]. Point vortices provide a fundamental example of systems with 
long-range interactions [8] which deserves a particular attention. We have shown that the main 
features of the kinetic theory: kinetic equation describing the evolution of the system as a 
whole, diffusion coefficient, drift term, Fokker-Planck equation describing the evolution of a 
test particle in a bath... could be obtained from a simpler formalism than the one developed 
in our previous papers [5U [23] . This clarifies the argumentation and delineates the domain of 
validity of the theory. We have given general equations that are valid for flows that are not 
necessarily axisymmetric nor markovian. A limitation of our approach is to neglect collective 
effects. These effects have been taken into account in [201 157] for axisymmetric flows. In plasma 
physics, collective effects are important because they lead to Debye shielding and regularize the 
logarithmic divergence at large scales that appears in the Landau equation (as shown by Lenard 
[21] and Balescu [22]). For point vortices, their influence seems less crucial since the kinetic 
equation (14"8|) derived by neglecting collective terms does not present any divergence. In addi- 
tion, concerning the expressions of the diffusion coefficient and drift term, we have indicated 
that collective effects have a negligible contribution because the diffusion coefficient and the 
drift velocity are dominated by local interactions. In future works, we plan to study in more 
detail the kinetic equations given in this paper. This project has been initiated in [29]. The 
kinetic theory could be used to interprete the numerical simulations of point vortices in 2D 
hydrodynamics [5S1 ES] or the experiments of non-neutral plasmas under a strong magnetic 
field (leading to quasi stationary states, vortex crystals,...) [50]. In agreement with the kinetic 
theory, these systems exhibit a violent collisionless relaxation and a slow collisional evolution. 
The collisionless relaxation is based on the 2D Euler equation and the evolution of the coarse- 
grained vorticity is described by Eq. (1841) . The collisional evolution was the main object of 
interest of the present paper. It is described by a general kinetic equation of the form (1521) that 
simplifies in Eq. ( I48p for axisymmetric flows. 
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A The potential of interaction 

The velocity of the i-th vortex is produced by the other vortices according to the relation 

Yt = --z x WiH = VV(j i), (149) 

where H is the Hamiltonian ([2]). The velocity created by point vortex j on point vortex i is 

V(j - <) = "7* x -p-. (150) 

Introducing a system of polar coordinates to localize the point vortices, the radial component 
in the direction of ri of the velocity created by point vortex 2 on point vortex 1 is 

VW2^1) = 2** (151) 
ri Wi 

In an infinite domain, the potential of interaction can be written 

«i2 = w(y rf + r%- 2r 1 r 2 cos(p) = u(r 1 ,r 2 , (f>), (152) 
where <fi = 6\ — 62. This implies that 

K 2 (1^2) = -^y n (2^1). (153) 

This relation results from the conservation of the angular momentum (see Appendix D of |23j) 
and remains valid in a bounded circular domain. Since the function u(r\, r 2 , 0) is periodic with 
period 2tt, it can be decomposed in Fourier series. Thus, 

M (n,r 2 ,0) = J2 eim ^m(ri,r 2 ), (154) 

m 

with 



u m {r u r 2 ) = — / cos(n<j))u(r 1 ,r 2 ,(j))d(j). (155) 
27r Jo 

Using the decomposition (j!54p . we find that Eq. (j!51j) can be rewritten in the form 

V ri {2 -+ 1) = i 7 iVmu ro (ri I r 2 )e im ( 91 -^. (156) 

m 

The usual potential of interaction between point vortices is solution of the Poisson equation 

Am = -5(r). (157) 

In an infinite domain, we have 

U12 = -— ln|ri -r 2 |. (158) 

The Fourier transform of u(<f)) can be easily obtained by using, e.g., the identities given in 
Appendix El of [23] . We find 

1 (r \ |m| 

u m (ri,r 2 ) = -—- ^ , (m^O) (159) 
47r|m| \r> / 
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Uo(ri,r 2 ) 



2vr 



lnr 



(160) 



where r> (resp. r<) is the largest (resp. smallest) of T\ and r 2 . In that case, the function 
defined in Eq. ( )47|) takes the explicit form 



x{n,r 2 ) 



—y 



1 / r< 



m=l 



< 



In 



(161) 



When the point vortices are confined within a circular box of radius R, the potential of inter- 
action is 



„ 12 = _J_ ln | ri _ r2 i + 



i? 2 
r x =-r 2 

To 



and its Fourier transform is given by 

1 



Ait m V r 



r> \ 2|m| 



(m ^ 0) 



(162) 



(163) 



I , ( R 2 
In 



27T \ r 2 r > 

In that case, the function defined in Eq. ( 1471) takes the explicit form 

2" 



x(n,^ 2 ) 



In 



+ 



4vr 2 



In 



8tt 2 



In 



(164) 



1 - 



r<r> 



(165) 



Finally, in the QG model, the potential of interaction is solution of an equation of the form 

1 



Am - -^-u = -S(r), 
L z 

where L is the Rossby radius. In an infinite domain, we obtain 



1 i<r (\ Yl ~ r2 l 
Ul2 = ^ K °{^- 



(166) 



(167) 



B A heuristic kinetic equation 



In [23] and in Sec. 13. 3[ we have derived a general kinetic equation (|52|) or (|54j) that is valid at 
order 0(l/iV). This equation can be simplified for axisymmetric flows (leading to Eq. (1481) ) 
and uni-directional flows (leading to Eq. (135) of [23]). In [23], we have heuristically proposed a 
generalized kinetic equation (137) that encompasses both the axisymmetric and unidirectional 
forms. This equation is not exact, so it cannot be obtained rigorously from Eq. fl52l) . Yet, it 
possesses interesting properties (conservation of E, T, L, P and if-theorem S > 0) so it can 
be useful. For axisymmetric and unidirectional flows, it does not exactly reduce to Eq. (T48l) 
and Eq. (135) of [23], but it has a similar structure so that the disagreement is not too severe. 
In this Appendix, we try to justify this equation but we stress that, since this equation is not 
exact, some approximations are necessarily un-controlled. 
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First, assuming that the decorrelation time is extremely short (which does not need to be 
the case) and that V ~ V, we replace Eq. ( |54|) by 

+ ^V>^ = 4 r* / W(2 M)V(2 -> 1,4 -t) 







^Jr^H 1 '^ 2 ' ' (168) 



where, now, the vorticity and vorticity gradient are evaluated at ri = ri(t) and r 2 = r 2 (t). 
Introducing the Fourier transform of the velocity created by point vortex 2 on point vortex 1: 

V(2 -)• l,t) = -ij J k ± n(A;)e ik(ri(t) - r2W) rfk, (169) 

and making a linear trajectory approximation Yi(t — r) ~ *i(t) — (V)(rj,t)r, we get 



d 3 



w( ri ,t)-(r 2 ,t),(170) 



<9r^ <9?"2 y ' 7 

where £ = r 2 — ri and v = (V) (r 2 , t) — (V) (ri, £) . The linear trajectory approximation is clearly 
not justified for an axisymmetric flow (since the point vortices follow circular trajectories as 
considered in Sec. 13 . 2j) so it again relies on the (un-controlled) hypothesis that the decorrelation 
time is extremely short. Integrating on time, we obtain 



^1 + ^(V}^ = -* 7 2 A J dvij <ik<ik'^*l"fi(< ; )fl(fc>-'< k+k ')^(k' ■ V) 



K^-£H ri ' i) > i) - (171) 

Now, using the heuristic argument that k' ■£ ~ 1, i.e. k' ~ £/£ 2 (due to the exponential e -ik '^), 
we make the rough substitution 5 (It' ■ v) — > A<5((£/£ 2 ) • v) = A£ 2 <5(£ ■ v), where A is a constant 
of order unity, in Eq. (11711) . Then, using Eq. (11691) . we can rewrite Eq. (|17ip in the form 

9Ul ' N ~\y)^-^ [dr 2 V»(2^ M)V(2-1, t)^(*.v) 



dt N N ' dr x dr{ 



x 



G^-a*)"^ (172) 



Finally, using V(2 — > 1) = — (7/27r)£j_/£ 2 , we obtain 



dt N 



This is Eq. (137) of [23]. The constant A has been determined so that Eq. (1 173ft reproduces 
at best the exact equations (Eq. ( )48i) and Eq. (135) of [23]) obtained for axisymmetric and 
unidirectional flows (in particular the expressions of the drift and diffusion in the corresponding 
Fokker-Planck equation, see [23]). This yields A = n/2 which is of order unity as expected. 
The previous arguments give some "justification" to Eq. (137) of [23] although we again stress 
that this equation is not exact so it is obtained from un-controlled approximations that are not 
really justified. 
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C Drift term for axisymmetric systems 

For an axisymmetric flow, the linearized equation (11 101) for the perturbation becomes 

d5u 1 d5ui T T , „ , oL> , 

■w^-rw = - K(F ^ 0) *- (174) 

Introducing the angular velocity and the potential of interaction, it can be rewritten 

d5uJ + 0( r ) — - 7 9MoP du f 1 7^ 

Taking the Fourier transform of this expression with respect to the angular variable, we obtain 

d8u} m . 7. * , \ do; 

— — + zmftto m = -^imu m (r,r P )e~ imdp — . 176 
at r ar 

Integrating this first order differential equation with respect to time, we find that 

Su m (t) = - f dre-* mnT ^mu m (r, r P ) e - ime ^)^. (177) 
J r ar 

To leading order in N — > +00, the test particle is advected by the mean flow so that 9p(t — r) = 
6p — QpT. Extending the time integral to +00, we obtain 

Su m (t) = - dTe- im ^- np ^tmu m (r J r P )e- im9p ^. (178) 
Jo r dr 

The integral on time r can be easily calculated yielding 

6u> m (t) = -n5\m(tt - ttp)]^imu m (r,r P )e- imep( ^. (179) 

r dr 

The drift velocity experienced by the test particle is 

ypol = i f ^ p^ffidTx, (180) 

7 J 



or, more explicitly, 



"ndn r ( w 1 —^6u 1 . (181) 



Introducing Fourier transforms and substituting Eq. ( 1179[> in Eq. (11811) . this can be rewritten 

KPO / = 2v^7 /' +0O dri y- m 2^ m(ri)rp) 2 5[m(ni _ np)] ^i (lg2) 
r P J dr x 



Finally, introducing the notation (147j) . we obtain 



= ^1 f + °° dr lX (r u r P )S(Q 1 - n P )^. (183) 
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